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W : ABSTRACT 

Nadarajah and Gupta (2004) introduced the beta Frechet (BF) distribution, which is a ge- 
neralization of the exponentiated Frechet (EF) and Frechet distributions, and obtained the 
probability density and cumulative distribution functions. However, they do not investigated 
its moments and the order statistics. In this paper the BF density function and the den- 
sity function of the order statistics are expressed as linear combinations of Frechet density 
^ ! functions. This is important to obtain some mathematical properties of the BF distribution 
I in terms of the corresponding properties of the Frechet distribution. We derive explicit ex- 
pansions for the ordinary moments and L-moments and obtain the order statistics and their 
moments. We also discuss maximum likelihood estimation and calculate the information 
^ . matrix which was not known. The information matrix is easily numerically determined. 

OO ! Two applications to real data sets are given to illustrate the potentiality of this distribution. 
O 

Keywords: Beta distribution. Exponentiated Frechet, Frechet distribution. Information 
KJi matrix, Maximum likelihood estimation. 

: 1. INTRODUCTION 

The Frechet distribution has applications ranging from accelerated life testing through to 
earthquakes, floods, horse racing, rainfall, queues in supermarkets, sea currents, wind speeds 
and track race records. Kotz and Nadarajah (2000) give some applications in their book. 
In this paper, we discuss the BF distribution which stems from the following idea. Eugene 
et al. (2002) defined the beta G distribution from a quite arbitrary cumulative distribution 
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function (cdf) G{x) by 



where a > and 6 > are two additional parameters whose role is to introduce skewness 
and to vary tail weight and B{a, b) = uj"'^^{l — uj^^^duo is the beta function. The class of 
distributions ([1]) has an increased attention after the works by Eugene et al. (2002) and Jones 
(2004). Application oi X = G~^{y) to the random variable V following a beta distribution 
with parameters a and 6, ~ -B(a, h) say, yields X with cdf ([T]). 

Eugene et al. (2002) defined the beta normal (BN) distribution by taking G{x) to be the 
cdf of the normal distribution and derived some of its first moments. General expressions for 
the moments of the BN distribution were derived (Gupta and Nadarajah, 2004). Nadarajah 
and Kotz (2004) also introduced the beta Gumbel (BG) distribution by taking G{x) to be the 
cdf of the Gumbel distribution and provided closed-form expressions for the moments, the 
asymptotic distribution of the extreme order statistics and discussed the maximum likelihood 
estimation procedure. Nadarajah and Gupta (2004) introduced the BP distribution by taking 
G{x) to be the Frechet distribution, derived the analytical shapes of the probability density 
function (pdf) and the hazard rate function and calculated the asymptotic distribution of 
the extreme order statistics. However, they do not investigate expressions for its moments 
and the information matrix which we do in this paper. Also, Nadarajah and Kotz (2005) 
worked with the beta exponential (BE) distribution and obtained the moment generating 
function, the first four cumulants, the asymptotic distribution of the extreme order statistics 
and discussed the maximum likelihood estimation. We can write ([1]) as 

F(x) = /G(.)(a,6), (2) 

where Iy{a,b) = B{a,b)^^ w"-^^{l — w)^^^dw denotes the incomplete beta function ratio, 
i.e., the cdf of the beta distribution with parameters a and h. For general a and 6, we can 
express (|2]) in terms of the well-known hypergeometric function defined by 

where (a)j = a{a + 1) . . . (a + i — 1) denotes the ascending factorial. We obtain 

a B[a, 0) 

The properties of the cdf F{x) for any beta G distribution defined from a parent G{x) in 
([1]), could, in principle, follow from the properties of the hypergeometric function which are 
well established in the literature; see, for example. Section 9.1 of Gradshteyn and Ryzhik 
(2000). 

The probability density function (pdf) corresponding to ([1]) can be written in the form 

m = ^^fi^r'i^ - Gix)}''-'gix), (3) 
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where g{x) = dG{x)/dx is the pdf of the parent distribution. The pdf f{x) will be most 
tractable when the functions G{x) and g{x) have simple analytic expressions as is the case 
of the Frechet distribution. Except for some special choices for G{x) in ([T]), it would appear 
that the pdf f{x) in ([3]) will be difficult to deal with. 

The cdf and pdf of the Frechet distribution are, respectively, 



a,A(x) = e-(^)\ x>0, (4) 

and 

g^,^(x) = Aa^x-(^+i)e-(^)\ a; > 0, 

where a > is the scale parameter and A > is the shape parameter. The rth moment of 
the Frechet distribution for r < X is fi[, = a'"r{l — r/X), and then the first four cumulants if 
A > 4 are 

^ 5^3 - 3^i5'2 + '^gl ^ Qa - 4^9193 + Q9I92 - ^gf 

i92-9!r^' ' (92-9!)' 

where gk = r(l — /c/A) for /c = 1, . . . , 4. 

Nadarajah and Gupta (2004) give the cdf of the BF distribution with parameters a > 0, 
6 > 0, cr > and A > (denoted by BF(a, b, a, A)) in the same way from ([1]) by replacing 
the parent cdf G{x) by (jl]) 

' .a-m , .\b~l. 



^^^) = 7?7-^/ ^"-^(l-^)''-M^ = /exp{-(^)A}(«,&), x>0. (5) 

li[a, 0) Jq 

They also give the corresponding pdf and hazard function, respectively, as 

fix) = -^^x-^'^'^e'^^^^'{l-e-&r~\ x>0, (6) 
B{a, h) 

and 

Aa^x-(^+i)e-«(^)'{l-e-(^)'}^-i 

= u( MJ1 — ( A ' > °- 

B{a, b){l - /exp{-(^)A}(a, b)} 

Figures [1] and [2] illustrate some of the possible shapes of the pdf and hazard function ([7]), 
respectively, for selected parameter values, including the case of the Frechet distribution. The 
BF distribution is easily simulated from ([5]) as follows: if V ~ B{a, b) then X = a/ (— log\^)^/^ 
has the BF(a, 6, cr. A) distribution. 

The BF distribution generalizes some well-known distributions. The exponentiated Frechet 
(EF) distribution (Nadarajah and Kotz, 2003) is a special case when a = 1. The Frechet 
distribution (with parameters a and A) is also a special case of ([6]) when a = b = 1. Further, 
when 6=1 and A = 1, ([6]) is an inverse gamma distribution with shape parameter 2 and 
scale parameter aa. Since the BF distribution generalizes the Frechet and EF distributions 
by adding two parameters and one parameter, respectively, it can be used by practitioners 
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Figure 1: Pdf of the BF distribution for selected parameter values. 

as an extra tool to analyze the data we would normally use with the last two distribu- 
tions. The book of Kotz and Nadarajah (2000) demonstrates the applicability of the Frechet 
distribution in several fields. 

The rest of the paper is organized as follows. Section 2 gives expansions for the pdf and 
cdf of the BF distribution and for the density of the order statistics depending on whether 
the parameters b (or a > 0) is real non-integer or integer. We show that the density functions 
of the BF and its order statistics can be expressed as mixture of Frechet density functions. 
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X=1 and 0=1 



X=1 and 0=2 



Fr6chet 
a.0.5,b.0.5 
a.0.5,b.2 
a.2,b.0.5 



Fr6chet 
a.0.5,b.0.5 
a.0.5,b.2 
a.2,b.0.5 



X=2 and 0=1 



X=2 and 0=2 



— Fr6chet 
a.0.5,b.0.5 
a.0.5,b.2 

- - a.2,b.0.5 



— Fr6chet 
a.0.5,b.0.5 
a.0.5,b.2 

- - a.2,b.0.5 



Figure 2: Hazard function of the BF distribution for selected parameter values. 



The moments of this distribution and of the order statistics are not known and general ex- 
pansions are derived in Section 3 for the cases 6 > real non-integer and integer. L-moments 
(Hosking, 1986) are expectations of certain linear combinations of order statistics and form 
the basis of a general theory which covers the summarization and description of theoretical 
probability distributions. In Section 4 we present expansions for the L-moments of the BF 
distribution. We discuss in Section 5 maximum likelihood estimation and calculate the ele- 
ments of the information matrix. Section 6 provides two applications to real data sets. In 
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Section 7 we end with some conclusions. 



2. EXPANSIONS FOR THE DISTRIBUTION AND DENSITY FUNCTIONS 

Here, we provide simple expansions for the cdf of the BF distribution depending on 
whether the parameter b (or a) is real non-integer or integer. We consider the series expansion 



valid for \z\ < 1 and b > real and non-integer. Application of ([8]) to ([T]) if 6 is real 
non-integer gives 

r(a) 2^Jib-j)j\{a + jy 
where Gcr,x{x) comes from (jlj). Then, we have 

r(. + i)^ (-iye-<-^^'(i)' 
rw ^r(6-j)j!(a + jr ' ' 

For 6 integer, the sum in (fTO!) simply stops at 6—1. When 6 = 1, it follows F{x) = e'""^^^^ . 
It can be seen in the Wolfram Functions Site @ that for integer b 

and for integer a 

Hence, if b is integer, we obtain another equivalent form for flTU]) 

,_e-(f)-^ r(a + j)(l-e-(f)^y 
r(a) ^ ,! 



and, for integer values of a, we have 



^http : //functions . wolf reun . com/ 
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If a = 1, the above expression reduces to 

F{x) = l-{l-e-(f)'}^ 

which agrees with the cdf of the EF distribution. 

The pdf in is straightforward to compute using any statistical software. However, 
we show that the BF density can be expressed as an infinite (or finite) weighted hnear 
combination of pdf 's of random variables having Frechet distributions. This is important to 
provide some mathematical properties of the BF distribution directly from the corresponding 
properties of the Frechet distribution. If 6 > is real non-integer, and again using ([H]) we 
can rewrite ([HD as 

oo 

f{x) =^Wkga,,x{^), (11) 

A;=0 

where 

Wk = r(a + 6)(-l)V{r(a)r(6 - k)k\{k + a)} 

represent weighted constants such that Ylh=o'^k ~ ^ 9ak,x{^) is a Frechet density with 
scale parameter = c^k + aY^'^ and shape parameter A. In addition, if a = 1, fill I) agrees 
with the corresponding result obtained by Nadarajah and Kotz (2003, Section 5). If 6 > 
is integer, the sum in ffTTl) is finite and stops at 6 — 1. Then, the ordinary, central, factorial 
moments and the moment generating function of the BF distribution could in principle 
follow from the same weighted infinite (or finite if b is integer) linear combination of the 
corresponding quantities for the Frechet distribution. 

We now give the density of the ith order statistic Xj:„, fi:n{x) say, in a random sample 
of size n from the BF distribution. It is well known that 

for i = 1, . . . ,n. Using ([5]) and ([6]) we can express fi:n{x) in terms of the incomplete beta 
function ratio 

The cdf of the ith order statistic Xi-_n, Fi,n{x) say, is 

Fi:n{,x) = ^ ( J/G;,,„(^)(a,6)''J{i_G,.„(^)}(6,a)""''. 

r=i ^ ^ 

Using the identity iYlT=o'^kx'^)"' ~ J2T=o^k,nx'' (see Gradshteyn and Ryzhik, 2000), 
where co,n = 

k 

Ck,n = {kaor^ - k + l)aiCk-i,n 

1=1 
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for A; = 1, 2, . . . and (fTOl) . the pdf of the ith order statistic can be written for 6 > real 
non-integer and integer, respectively, as 

, , ,_^^ i-^nV)nby^'-'Biai^ + k)+J,b)c^l ^ 

and 

" * i~^)''Ck^)B{aii + k) + j,b)c^-^\ 

fi:n{x) = 2, y, TT? u\j+k TDl ■ • , 1 \ '~fi,j,k{x), (13) 

where fij^k{x) is the denstity of a BF(a(2 + /;;)+ j, b, a, A) distribution. 



(1) ^J^V^'"' _^r(&).A (-iy{/(. + A:)-j} 



and 

(2) _/iy+'=-^ _a^ (-iy(V){/(^ + fc)-j} 

for j > 1. Expansions (fT0|) - (fT3|) are the main results of this section. 

Two alternative expansions for the densities of the order statistics follow from the identity 
(Zli^i (^i)^ = Z]{^i,...,mfc}=o «mi • • • CLnik ^OT k a positivc integer. Using this identity and ([ID]), 
it is easy to show for 6 > real non-integer and integer that 

n—i oo oo 

/-(^) = EE--- E ^Sm^) (14) 

k=0 mi=0 mi_|.fe_i=0 

and 

n—i b—1 b—1 
A;=0 mi=0 mi_|_fe_i=0 

respectively, where fi^k{x) is the pdf of a BF(a(i + k) + X]j=i"^ "^j; ^) distribution, 

^(1) ^ (-i)^+gr^i-^ (T)^(»(^ + ^) + b)T{br^~' 

"'^ B{a, bY+^B{i, n-i + l) IljS"^ r(6 - mj)mj!(a + m^) 

and 

_ (-1)^+^"-'^ ("r)i^(a(^ + fe) + Egr^ b) '1^' (t;) 



B{a,by+''B{i,n-i + l) (a + mj)' 



The summation in (fT^ and (fT5|) extends over all {i + A;)-tuples {k,mi, . . . ,mj_|_fc_i) of non- 
negative integers and is implementable on a computer. However, expansions (1121) and fll3l) 
are much simpler to be calculated and their CPU times are usually smaller. 



3. MOMENTS 

As with any other distribution, many of the interesting characteristics and features of 
the BF distribution can be studied through the moments. We obtain immediately the rth 
moment /i^ of the BF distribution from (ITT!) if r < A 

= rw g ra-m ■ <«) 

If 6 > is integer and r < A, the sum stops at 6 — 1. If a = 1 and r < A, (fT6l) gives the rth 
moment of the EF distribution with parameters b, a and A which is a new result for the EF 
distribution. 



From (1121) and (jT3l) . we obtain simple expansions for the moments of the order statistics. 
The rth moment of the for 6 > real non-integer is 

:^^(-l)H"7')T(by+''-^B(a(t + k)+ jMc^P,^ 
V i.n/ Bia,hY+^Bii,n-i + l) ^ ^ ' 

k=0 j=0 \ ^ J V ' 'J 



and for 6 > integer 

n—i b—1 / 1 

B{a,by+''B{i,n-i + l) 



fc=o i=o 



where Xj^^fc ~ BF(a(i + k) + j,b,a, X). From (HM and (fTSl) . we obtain two alternative 
expressions for the moments of the order statistics. These expressions for 6 > real non- 
integer and integer are given by 



n—i oo 



k=0 mi=0 mi_|_fe_i=0 



and 

n—i b—1 b—1 



k=0 mi=0 ?Tii+fe_i=0 



respectively, where Xi^k ~ BF{a{i + k) + X]!=i ^ "^i^ ^5 ^! •^)- 
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1.3 
b 



Figure 3: Skewness of the BF distribution as a function of a{b) for selected values of b{a). 




Figure 4: Kurtosis of the BF distribution as a function of a{b) for selected values of b{a). 

Graphical representation of skewness and kurtosis when A = 5 and o" = 1, as a function 
of parameter a for some choices of parameter b, and as a function of parameter b for some 
choices of parameter a, are given in Figures 3 and 4, respectively. 
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4. L-MOMENTS 



The L-moments are analogous to the ordinary moments but can be estimated by hnear 
combinations of order statistics. They are hnear functions of expected order statistics defined 
by (Hosking, 1990) 



Xr 



(r + 11 



k=0 



(19) 



The first four L-moments are: Ai = E{Xi.,i), A2 = lE{X2:2 - Xi,2), A3 = lE{Xs:3 - 2X2^3 + 
Xi-s) and A4 = jE{X4,4 — 8X3.4 + 8X2:4 — X1.4). The L-moments have the advantage that 
they exist whenever the mean of the distribution exists, even though some higher moments 
may not exist, and are relatively robust to the effects of outliers. 

From the expansions (fT7|) - (fT8|) for the moments of the order statistics we can obtain ex- 
pansions for the L-moments of the BF distribution as weighted linear combinations of the 
means of suitable BF distributions. 



5. ESTIMATION AND INFORMATION MATRIX 

We assume that Y follows the BF distribution and let 6 = (a, b, a, X)^ be the true 
parameter vector. The log-likelihood i = i{6) for a single observation y of Y is given by 



log A + A log{a/y) - log{5(a, b)} - a{a/yf + (6 - 1) log{l - e 



}• 



The components of the score vector U 
observation are given by 

da 
d£ 

db 

m 
m 

dX 



U{e) = {di/da, di/db, di/da, di/dXf for one 
-V'(a) + ^(a + 6) - {(T/y)\ 



-ipib) +^(a + 6) + log{l 
A A(T^-i r 6-1 



-i-Y 
e y 



y 



{S-)\ 

e y 



Y + log(- 

A \y 



1 



b-l 

e y — 



From E{di/da) = 0, we obtain 



E{X 



-An ip{a + b) - ip{a) 



For interval estimation and hypothesis tests on the model parameters, we require the infor- 
mation matrix. The 4x4 unit information matrix K = K{9) is 

^ ^a,a ^a,6 /^a.o- '^a,A ^ 



K 



l^a,b l^b,b l^b.cr l^b,\ 
^a,a ^b.a ^a,a ^a,X 
\ l^a,X l^b,X l^cr,X f^X,X ) 
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whose elements are 



i^a,a = ^[l + a(A-l){V'(a + 6)-V'(a)} + (6-l)(Ari,i,2,o-Ti,i,i,o)], 

«A,A = ^{1 + 0^0,0,1,2 + (^'-l)(2'l,2,2,2-Ti,i,i,2)}, 

i^a,x = --[1 -a{V'(a + 6) -V'(a)+ro,o,i,i} + (^'-l)(T^i,i,i,o + 7^1,1,1,1 

(T 

-ATi,2,2,o)], 'ta.b = -Ip'ia + 6), Ka^A = ^^0,0,1,1' 
«a,<T = -{^(a + - V'(a)}, K6,a = --71,1,1,0, «6,A = -^^1,1,1,1- 

a a A 

Here, we define a random variable V following a Beta{a, b) distribution and the expected 
value 

= E [v\i - y)-^(-iogy)Hiog(-iogy)}'] , 

where the integral obtained from the above definition is numerically determined using MAPLE 
and MATHEMATICA for any a and h. For example, for a = 1.5 and b = 2.5 we easily 
calculated all T's in the information matrix: Ti,i,2,o = 0.51230070, Ti,i,i,o = 0.55296103, 
7o,o,i,2 = 0.62931802, Ti,2,2,2 = 0.43145336, Ti,i,i,2 = 0.32124774, To,o,'i,'i = 0.48641180, 
Ti'i'i'i = -0.16152763 and Ti,2,2,o = 0.86196008. ' " 

For a random sample y = {yi, . . . , of size n from Y, the total log-likelihood is 



i=l 



where fS^^ is the log-likelihood for the ith observation {i = 1, . . . ,n) as given before. The 
total score function is Un = Un{9) = Y17=i U^'^\ where U^'''^ for i = 1, . . . ,n has the form 
given earlier and the total information matrix is Kn{6) — nK{d). 

The MLE ^ of ^ is numerically determined from the solution of the nonlinear system of 
equations f/„ = 0. Under conditions that are fulfilled for the parameter 9 in the interior 
of the parameter space but not on the boundary, the asymptotic distribution of -Jn{9 — 
9) is NiiQ, K{9y^). The asymptotic multivariate normal N^{Q, Kn{9)~^) distribution of ^ 
can be used to construct approximate confidence regions for some parameters and for the 
hazard and survival functions. In fact, an 100(1 — 7)% asymptotic confidence interval for 
each parameter 9i is given by 

ACh = {9i - z^/2^k^% 9i + z^/2Vk^), 

where k,^^'^* denotes the ith diagonal element of Kn{9)~^ for i = 1,2,3,4 and is the 
quantile 1 — 7/2 of the standard normal distribution. The asymptotic normality is also 
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useful for testing goodness of fit of the four parameter BF distribution and for comparing 
this distribution with some of its special submodels using the likelihood ratio (LR) statistic. 

We consider the partition 6 = (9j, 9j)^ , where 6i is a subset of the parameters of interest 
of the BF and 62 is a subset of the remaining parameters. The LR statistic for testing the 
null hypothesis Hq : di — df^ versus the alternative hypothesis Hi : di ^ Of^ is given by 
w — 2{^(^) —l{9)}, where 9 and 9 denote the MLEs under the null and alternative hypothe- 
ses, respectively. The statistic w is asymptotically (as n 00) distributed as x^,, where 
k is the dimension of the subset 9i of interest. Then, we can compare for example a BF 
model against an EF model by testing Hq : a = 1 versus Hi : a ^ 1. We can also com- 
pare a BF model against the Frechet model by testing Hq : a — h — 1 versus Hi : Hq is false. 

6. APPLICATIONS 

In this section we fit the BF distribution to two examples of real data and test two types 
of hypotheses: Hq : Frechet x Hi : BF and Hq : EF x Hi : BF. The first example is 
an uncensored data set from Nichols and Padgett (2006) consisting of 100 observations on 
breaking stress of carbon fibres (in Gba): 3.7, 2.74, 2.73, 2.5, 3.6, 3.11, 3.27, 2.87, 1.47, 3.11, 
4.42, 2.41, 3.19, 3.22, 1.69, 3.28, 3.09, 1.87, 3.15, 4.9, 3.75, 2.43, 2.95, 2.97, 3.39, 2.96, 2.53, 
2.67, 2.93, 3.22, 3.39, 2.81, 4.2, 3.33, 2.55, 3.31, 3.31, 2.85, 2.56, 3.56, 3.15, 2.35, 2.55, 2.59, 
2.38, 2.81, 2.77, 2.17, 2.83, 1.92, 1.41, 3.68, 2.97, 1.36, 0.98, 2.76, 4.91, 3.68, 1.84, 1.59, 3.19, 
1.57, 0.81, 5.56, 1.73, 1.59, 2, 1.22, 1.12, 1.71, 2.17, 1.17, 5.08, 2.48, 1.18, 3.51, 2.17, 1.69, 
1.25, 4.38, 1.84, 0.39, 3.68, 2.48, 0.85, 1.61, 2.79, 4.7, 2.03, 1.8, 1.57, 1.08, 2.03, 1.61, 2.12, 
1.89, 2.88, 2.82, 2.05, 3.65. 

The MLEs and the maximized log-likelihood using the BF distribution are 
d = 0.4108, b = 125.1891, A = 0.7496, a = 31.4556, Ibf = -142.9640, 
whereas for the EF and Frechet distributions we obtain 

b = 52.0491, A = 0.6181, a = 26.1730, Ief = -145.0870, 

and 

A = 1.7690, a = 1.8916, %^chet = -173.1440, 

respectively. 

The second data set is obtained from Smith and Naylor (1987). The data are the strengths 
of 1.5 cm glass fibres, measured at the National Physical Laboratory, England. Unfortu- 
nately, the units of measurement are not given in the paper. The data set is: 0.55, 0.93, 
1.25, 1.36, 1.49, 1.52, 1.58, 1.61, 1.64, 1.68, 1.73, 1.81, 2 ,0.74, 1.04, 1.27, 1.39, 1.49, 1.53, 
1.59, 1.61, 1.66, 1.68, 1.76, 1.82, 2.01, 0.77, 1.11, 1.28, 1.42, 1.5, 1.54, 1.6, 1.62, 1.66, 1.69, 
1.76, 1.84, 2.24, 0.81, 1.13, 1.29, 1.48, 1.5, 1.55, 1.61, 1.62, 1.66, 1.7, 1.77, 1.84, 0.84, 1.24, 
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1.3, 1.48, 1.51, 1.55, 1.61, 1.63, 1.67, 1.7, 1.78, 1.89. 



Fitting the BF, EF and Frechet distributions we obtain the MLEs and the maximized 
log-hkehhood: 

a = 0.3962, b = 225.7272, A = 6.8631, a = 1.3021, Ibf = -90.5180, 
b = 112.5986, A = 7.7859, a = 0.9814, Ief = -93.1962 

and 

A = 1.2643, a = 2.8875, %echet = -117.7765, 

respectively. 

For the first data set, the values of the LR statistics for testing the hypotheses Hq : 
Frechet x Hi : BF and Hq : EF x Hi : BF are: 60.36 (p-value=7.81 x 10~^^) and 4.246 
(p-value=3.93 x 10~^), respectively. For the second data set, we obtain the values of the LR 
statistics 54.5170 (p-value=1.45 x 10~^^) and 5.3564 (p-value=2.06 x 10~^) for the hypotheses 
Hq : Frechet x Hi : BF and Hq : EF x Hi : BF, respectively. Therefore, in both situations, 
using any usual significance level we reject the null hypotheses in favor of the alternative 
hypothesis that the BF distribuiton is an adequate model. 

The plots of the estimated densities of the BF, EF and Frechet distributions given in 
Figure [5] show that the BF distribution gives a better fit than the other two submodels for 
both data sets. 

7. CONCLUSIONS 

The BF distribution provides a rather general and flexible framework for statistical ana- 
lysis of positive data. It unifies some previously proposed distributions, therefore yielding a 
general overview of these distributions for theoretical studies, and it also provides a rather 
flexible mechanism for fitting a wide spectrum of real world data sets. The BF distribution 
is motivated by the wide use of the Frechet distribution in practice, and also for the fact that 
the generalization provides more flexibility to analyze more complex situations. In fact, the 
BF distribution ([6]) represents a generalization of some distributions previously considered 
in the literature such as the Frechet and EF (Nadarajah and Kotz, 2003) distributions. 
This generalization provides a continuous crossover towards cases with different shapes (e.g. 
skewness and kurtosis). 

The BF density can be expressed in the mixture form of Frechet densities. For doing 
this, we derived some expansions for the cdf of the BF distribution and their ordinary and 
L-moments. We call attention for the fact that the moments of the EF are not known in the 
literature and we derived these moments as a particular case of our results. The pdf of the 
BF order statistics can also be expressed in terms of a linear combination of Frechet densi- 
ties. We also derive the moments of the order statistics. We discuss the maximum likelihood 
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Figure 5: Fitted densities of the BF, EF and Frechet distributions for the data sets 1 and 2, 
respectively. 

estimation and obtain the information matrix, and considered the LR test which may be very 
useful in practice. We show that the formulae related with the BF are manageable, and with 
the use of modern computer resources with analytic and numerical capabilities, may turn 
into adequate tools comprising the arsenal of applied statisticians. Two numerical examples 
illustrate that the BF distribution provides better fits than the EF and Frechet distributions. 
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